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The SU(3) running coupling from lattice gluons 

UKQCD Collaboration - presented by C. Parrinello* 

Department of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3JZ, UK. 

We provide numerical results for the running coupling in 5(7(3) Yang-Mills theory as determined from an 
analysis of lattice two and three-point gluon correlation functions. The coupling is evaluated directly, from first 
principles, by defining suitable renormalisation constants from the lattice triple gluon vertex and gluon propagator. 
For momenta larger than 2 GeV, the coupling is found to run according to the 2-loop asymptotic formula. The 
influence of lattice artifacts on the results appears negligible within the precision of our measurements, although 
further work on this point is in progress. 



1. INTRODUCTION 

The running coupling as{q), where g is a mo- 
mentum scale, is a fundamental QCD quantity 
providing the link between the low and high- 
energy properties of the theory. Given a renor- 
malisation scheme, as{q) can be measured exper- 
imentally for a wide range of momenta. A pre- 
cise determination of as(g) (or equivalently of the 
scale A determining the rate at which a^ runs) is 
extremely important as it would fix the values of 
various parameters in the Standard Model, pro- 
viding bounds on new physics. 

Computing the QCD running coupling is a ma- 
jor challenge for the lattice community. In this re- 
spect, several different methods have so far been 
investigated [1-4]. Such methods, excepting stud- 
ies of the static potential, share the disadvan- 
tage that for each choice of the simulation pa- 
rameters only one value of the running coupling 
can be determined. Thus, in general, several dif- 
ferent data sets are needed in order to observe 
the scale dependence of the coupling. On the 
other hand, when extracting the coupling from 
the static potential, one has to fit the data to 
some phenomenological expression, so that the 
coupling is determined in a rather indirect way. 

Here we investigate the feasibility of a more 
recent proposal [5]. It differs from the previous 
ones in that it adopts a definition of the cou- 
pling directly from the fundamental fields appear- 
ing in the continuum Lagrangian. Such a method 



allows the determination of the coupling at dif- 
ferent scales from a single simulation, and also 
has the advantage that the extension to full QCD 
presents, in principle, no extra problems. 

2. THE METHOD 

Consider the standard lattice definition of the 
gluon fields in terms of the link variables [6], 
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traceless 

where a is the lattice spacing and go the bare 
coupling. One can define the n-point gluon cor- 
relation functions in momentum space as 

Gtl.APl^- • -.Pn) = (^Mi(Pl) • ■■A.APn)), (2) 

imposing momentum conservation. In particular, 

G%\p)=T,,{p)G{p') (3) 

is the standard gluon propagator. The amputated 
three-gluon vertex is defined as 
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G(2)(pi)G(2)(p2)G(2)(p3) ^^^ 

The vertex renormalisation constant, Zv{qa), 
can be defined at any fixed momentum scale q by 
choosing appropriate kinematics to project out 
the part of F proportional to the tree-level ver- 
tex [5]. We work in the Landau gauge, and not- 
ing the general form of F in the continuum (in 
any covariant gauge [7]), we evaluate F^i at the 
kinematical points defined by 



0^ = 1 i^P, Pl = -P3=Pp, P2 = 0. 



(5) 



Here, pp = q ep, where ep is the unit vector in 
the /? direction. Then, in the limit a —f 0, T can 
be written as 



rL/3a(-P^.0,p^) = 2 Zy\qa) go q. 



(6) 



Analogously, the wavefunction renormalisation 
constant Za{<1<i) is defined from the gluon prop- 
agator by 

G{q') = ZA{qa)^. (7) 

The renormalised coupling g[q) (where as{q) = 
£r^(g)/47r) is finally obtained via 

73/2 
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From our definition, one can check explicitly that 
g[q) is independent oi go (as it should be). 

3. COMPUTATIONAL PROCEDURE 

SU{S) gauge configurations were produced us- 
ing a hybrid-overrelaxed algorithm, every sixth 
sweep being a pseudo-heat-bath. We generated 
three data sets, each comprising 150 configura- 
tions, on a IS"* lattice. The values of /? = 6/gl 
were chosen as 6.0, 6.2 and 6.4. Successive con- 
figurations were separated by 150 sweeps, 1000 
sweeps being allowed for thermalisation. 

A crucial step in the method is the accurate 
implementation of the lattice Landau gauge con- 
dition 



A(a;) = ^^^(a; + e^)-^^(a;) = 0, Va 



(9) 



This is achieved by using a Fourier-accelerated 
algorithm [8]. To monitor the gauge-fixing accu- 
racy we compute the quantity 
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as the algorithm progresses, terminating when 
< lO""""""" (this is close to 32-bit machine pre- 
cision). Since our calculations involve low mo- 
mentum modes of gluon correlation functions, the 
above test is not in itself sufficient since A(a;) is a 
local quantity. For this reason, we also compute 



Ao{t) = J2M^,i)- 



(11) 



In a periodic box, the Landau gauge condition 
implies that Ao(t) is independent oit [6]. For our 
configurations, Ao(t) is constant to better than 
one part in 10^. 

For the purpose of our calculation, the only 
quantity which we need to store is the Fourier 
transform of j4^(a;) for a selected range of lattice 
momenta. All n-point gluon correlation functions 
can then be assembled at the analysis stage using 
Eq. 2, where momentum conservation is explic- 
itly imposed. The calculation of the renormalised 
running coupling then proceeds as follows. First, 
we evaluate the gluon propagator and compute 
Za from Eq. (7). Next, we compute the com- 
plete three-point function G^^' of the gluon field. 
We obtain F, the amputated vertex function, by 
dividing G^^) by the propagators corresponding 
to the momenta on the external legs. This yields 
Zv from Eq. (6). Finally, g[q) is obtained from 
Eq. (8). We fully exploit the symmetries of Eq. 
(6) to increase statistics. 

4. RESULTS 

We plot the running coupling versus momen- 
tum in Fig. 1. We express q in physical units 
using values for the inverse lattice spacing taken 
from string tension measurements [2,9]. 

We obtain a very clear signal despite the com- 
plicated numerical procedure. The important 
physical question is whether there exists a range 
of momenta for which the coupling runs according 
to the 2-loop perturbative expression 
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Here, 60 = ll/167r2, bi = 102/(167r2)2 and A 
is the QCD scale parameter for the pure gauge 
theory. To answer this question, and to obtain 
an estimate for A, we compute A as a function 
of the measured values of g[q) according to the 
formula 
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O /S = 6.0, a"'=1.94 GeV 
□ |S = 6.2, a"' = 3.73 GeV 
X ^ = 6.4, a"' = 3.62 GeV 
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Figure 1. Running coupling vs. momentum. 



Figure 2. Lambda parameter vs. momentum. 



If the coupling is running according to Eq. (12), 
then we expect A defined from the above equation 
to be constant. 

We plot A versus q in Fig. 2, and see that for 
g > 2 GeV the data are consistent with a constant 
value for the QCD scale parameter. 

Since we have made measurements on different 
physical volumes and with different lattice mo- 
menta qa, we are able to check that, within errors, 
our results are independent of both the infra-red 
and ultra-violet cutoffs. We only show results ob- 
tained with at most three units of lattice momen- 
tum. As expected, when using higher momenta 
(ga > 37r/8), finite lattice spacing effects become 
clearly visible. 

5. CONCLUSIONS 

We have shown that a non-perturbative deter- 
mination of the QCD running coupling can be 
obtained from first principles by a lattice study 
of the triple gluon vertex. Eventually, our values 
for the coupling g[q) (or equivalently our value 
for A) will be converted to predictions for the 
corresponding quantities in the MS scheme via 
a matching calculation. Also, we plan to further 
analyse the role of lattice artifacts. 
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